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Department  of  Chemical  Engineering  and  Fuel  Technology, 
Sheffield  University,  SI  3JD. 


Abstract 


Infrared  technology  finds  many  applications  in  the  field  of  combustion,  ranging  from  pollution  monitoring, 
through  military  systems,  to  the  control  of  industrial  furnaces  and  boilers.  This  review  of  some  selected 
concepts  highlights  the  interaction  between  the  diagnostic  role  of  infrared  measurements  and  the  current 
status  of  mathematical  modelling  of  combustion  systems.  The  link  between  measurement  and  computing  has  also 
evolved  to  the  point  where  a digital  processor  is  becoming  an  inherent  part  of  many  new  instruments.  This 
point  is  illustrated  by  reference  to  the  diffraction  particle  size  meter,  fire  detection  and  alarm  systems, 
and  furnace  control.  In  the  future,  as  fuels  become  scarce  and  expensive,  and  micro-electronics  become  more 
available  and  inexpensive,  it  is  certain  that  infra-red  devices  will  find  increasing  application  in  smaller 
industries  and  the  home. 


Int  roduction 


The  continuing  availability  of  energy  is  one  of  the  most  important  challenges  facing  mankind  at  the 
present  time.  Research,  development  and  demonstration  of  systems  which  avoid  waste  during  energy  conversion 
is  vital  if  we  are  to  successfully  meet  this  challenge.  By  the  year  2000,  it  has  been  estimated  that  con- 
servation measures  could  save  up  to  about  30Z  of  the  potential  national  fuel  requirements  with  no  loss  of 
'amenity  (D.  Although  nuclear  energy  will  contribute  approximately  a further  30Z , it  is  clear  that  com- 
bustion of  fossil  fuels  will  continue  to  represent  a major  part  of  the  energy  source.  In  the  short  term 

10  years)  it  is  apparent  that  oil  will  be  a major  fuel,  and  indeed  at  the  present  time,  about  half  the 
world’s  energy  is  derived  from  oil.  However,  the  oil  will  be  largely  exhausted  by  the  end  of  the  century 
whereas  the  coal  reserves  will  last  for  many  years.  It  therefore  follows  that  the  emphasis  in  fuel  tech- 
nology will  gradually  move  from  oil  to  coal  and  coal  derived  fuels. 

It  is  the  responsibility  of  scientists  and  engineers  to  undertake  the  R,D  and  D programs  which  will 
achieve  the  rational  use  of  energy,  and  instrumentation  is  the  key  to  all  phases  of  this  program.  Industrial 
applications  of  infrared  technology  in  the  energy  field  are  therefore  largely  concerned  with  combustion  and 
associated  conservation  techniques. 

In  this  paper,  it  is  not  possible  to  consider  all  the  relevant  applications  of  infrared  technology,  and 
therefore  only  a few  selected  concepts  will  be  reviewed.  These  items  have  been  chosen  in  the  hope  that  some 
at  least  will  be  novel  and  will  stimulate  further  work  in  this  important  and  exciting  field. 

It  is  currently  impossible  to  discuss  combustion  and  energy  without  including  the  minimization  of  pollu- 
tants as  an  integral  part  of  the  study.  It  is  noteworthy  that  chemical  pollutants  (CO,  N0X,  SOj , SO3,  dust, 
smoke)  are  two  or  three  orders  of  magnitude  more  dilute  than  the  species  of  interest  in  determining  com- 
bustion efficiency  (CO2,  O2)  and  instrumentation  systems  must  be  correspondingly  more  sensitive.  However, 
in  research  and  development,  it  is  not  usually  sufficient  to  simply  measure  the  properties  of  combustion 
systems,  but  rather  we  must  compare  the  results  with  the  predictions  of  levels  and  trends,  and  study  how 
these  vary  with  the  geometrical  and  operating  parameters  at  the  control  of  the  designer.  The  overall  result 
of  complementary  measurement  and  prediction  i6  the  rational  design  of  optimized  combustion  systems. 

In  addition  to  the  stimulus  to  infrared  measuring  techniques  in  combustion  caused  by  concern  over  energy 
and  pollution,  dramatic  innovations  have  been  brought  about  by  military  incentives.  The  heat  seeking  anti- 
aircraft missile,  infrared  signature  techniques  and  certain  aspects  of  earth  surveillance  sattelites  are 
obvious  examples  of  infrared  applications  involving  combustion. 

As  the  gradual  oil  shortage  forces  us  to  move  towards  fuels  with  a wider  specification  and  particularly 
coal  derived  fuels,  it  is  inevitable  that  the  hydrogen/carbon  ratio  of  the  liquid  fuel  will  decrease  from 
the  present  value  of  2 ('i  20Z  aromatics  for  aircraft  fuel)  to  a value  approaching  1.5.  This  change  will  tend 
to  cause  soot  to  form,  which  not  only  increases  the  smoke  pollution  problem  but  also  increases  radiant  heat 
transfer  leading  to  increased  loading  in  critical  boiler  regions  and  increased  flame  tube  temperatures  in 
gas  turbine  combustors. 

Combustcr  Modelling 


Presently  the  aviation  industry  produces  a wide  variety  of  engine  designs,  ranging  from  simple  lift 
engines  to  sophisticated  multispool  by-pass  engines.  In  so  doing  the  combustion  designer  has  the  important 
task  of  selecting  the  operating  conditions  which  yield  maximum  efficiency,  while  minimizing  the  emission  of 
pollutants.  Due  to  the  complexity  of  the  processes  involved  and  consequently  lack  of  fundamental  under- 
standing, the  designer  frequently  has  to  resort  to  general  rules-of-thumb , based  on  his  past  experience  and 
qualitative  analysis.  The  rationale  behind  mathematical  modelling  is  to  replace  this  rule-of-thumb  approach 
by  a firmer,  quantitative  procedure  while  optimizing  a combustor  design  for  any  given  application. 

The  technique  of  mathematical  modelling  has  been  the  outcome  of  two  disciplines  that  are  apparently  inde- 
pendent but  complementary  in  nature.  The  first  of  these  attempts  to  solve  the  differential  transport 
equations  for  the  reacting  turbulent  continuum  using  either  a finite  difference  or  a finite  element  tech- 
nique. This  lint  of  approach,  with  the  inclusion  of  proper  physical  models  and  powerful  computers,  provides 


superior  combustor  predictions  of  all  kinds.  The  accurate  local  prediction  of  two-phase  three-dimensional 
aerodynamic  effects  and  combustion  influences  is  now  feasible  via  such  models. 

Contrary  to  the  rigorous  description  employed  in  the  finite  difference  modelling  by  partial  differential 
equations  of  heat,  mass  and  momentum  transfer,  the  "Chemical  Reactor  Modelling"  approach  concentrates  on  the 
representation  of  the  combustor  flow  field  in  terms  of  interconnected  partially  stirred  and  plug  flow 
reactors.  This  approach  has  the  advantage  that  the  complex,  time-consuming  solution  of  the  equations  is  re- 
placed by  simple  flow  models  and  the  computational  requirements  are  generally  minimum.  These  features  enable 
the  attractiveness  of  chemical  reactor  modelling  to  be  exploited  during  the  development  phase  of  combustor 
design. 

The  problem  of  mathematical  modelling  combustors  is  complicated  by  the  fact  that  the  simultaneous  processes 
involved  comprise  three-dimensional  two-phase  fluid  dynamics,  turbulent  mixing,  fuel  evaporation,  radiative 
and  convective  heat  transfer  and  chemical  kinetics.  For  a model  to  claim  success  at  all  in  the  prediction  of 
such  inter-related  phenomena,  the  above  representative  processes  have  to  be  physically  modelled  using  the 
findings  of  previous  experimental  work.  The  problems  encountered  at  this  stage  of  modelling  stem  mainly  from 
the  complexity  of  the  laws  governing  turbulent  transport,  the  chemical  kinetics  of  laminar  and  turbulent  re- 
acting gases,  formation  and  disappearance  of  condensed-phase  particles  and  thermal  radiation  through  absorbing 
and  scattering  media.  The  other  source  of  difficulty  is  mathematical  in  nature  and  is  attributable  to  the 
fact  that  most  combustion  chamber  flows  are  three-dimensional  and  the  differential  equations  describing  them 
are  numerous,  simultaneous  and  non-linear. 

The  specific  combustor  modelling  strategy  adopted  in  Sheffield  comprises  a hybrid  approach  the  first 
stage  of  which  concentrates  on  the  prediction  of  hydrodynamic  and  thermal  patterns  via  a basic  finite  diff- 
erence scheme. ^ The  pertinent  features  of  turbulence  are  herein  simulated  by  models  incorporating  the  dis- 
tribution of  kinetic  energy  of  turbulence  and  its  dissipation  rate.  Spray  combustion  is  computed  from  the 
droplet  trajectories  and  evaporation  to  yield  concentrations  of  the  discrete  size  groups  of  droplets.  Chemical 
reaction  is  assumed  to  proceed  in  two  distinct  steps,  thereby  enabling  local  concentrations  of  CO,  C02,  Oo, 

*>2 » ^2®  and  fuel  to  be  evaluated. 

Due  to  the  attractiveness  of  the  flux  methods  in  combined  fluid  flow  and  heat  transfer  predictions,  radia- 
tion effects  (largely  infrared)  are  modelled  by  reference  to  the  six-flux  model  of  radiation.  The  different- 
ial equations  describing  the  variations  of  the  fluxes  are:- 
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ln  the  above  equations  the  various  terms  are  defined  as: 


cylindrical  polar  co-ordinates, 
flux  in  the  direction  of  positive  r, 

flux  in  the  direction  of  negative  r. 

flux  in  the  direction  of  positive  x. 

flux  in  the  direction  of  negative  x. 

flux  in  the  direction  of  positive  6. 

flux  in  the  direction  of  negative  6. 

flux  absorption  coefficient,  fraction  of  radiation  absorbed  per  unit  length, 
scattering  coefficient,  fraction  of  radiation  scattered  per  unit  length. 
oT**t  black  body  emissive  power  at  the  fluid  temperature. 

Stef an-Boltzmann  constant. 


Following  previous  studies,^  * the  numerical  solution  scheme  has  the  pressure  and  velocities  as  the  main 
flow  variables,  and  is  based  on  an  implicit,  upwind  finite  differencing  scheme.  (SIMPLE  algorithm).  The 
governing  equations  for  the  dependent  variables  are  integrated  over  finite  control  volumes  to  derive  the 
finite  difference  equations  which  are  then  solved  by  a procedure  proven  to  be  numerically  stable. 

Upon  obtaining  a compatible  set  of  hydrodynamic  and  thermal  fields,  the  second  phase  of  the  prediction 
procedure,  i.e.  chemical  reactor  modelling  is  initiated.  Here  attention  is  focussed  or,  combustion  efficiency 
and  pollution  levels  produced  by  a particular  combustor  design  using  an  extensive  17  step  kinetic  scheme.  The 
basic  criterion  employed  in  relating  the  two  stages  of  the  computations  is  to  recognize  that  mixing  at  the 
molecular  level  is  carried  out  by  the  movement  of  molecules  between  turbulent  eddies  of  different  concen- 
trations. Since  the  molecules  also  carry  momentum  and  hence  dissipate  turbulence,  it  is  hypothesized  that 
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mixing  rate  is  proportional  to  turbulence  dissipation  rate.  The  finite  difference  solution  yields  the  dis- 
tribution of  turbulence  dissipation,  high  values  prevailing  in  the  vicinity  of  the  gas  jet  injections.  These 
regions  thus  become  well  stirred  reactors  of  corresponding  volume  with  the  degree  of  stirring  or  mixing  rate 
factor  being  determined  from  the  total  dissipation  within  the  reactor.  Reactors  having  low  mixing  rates  are 
represented  as  plug  flow  reactors  and  it  is  assumed  that  reactions  taking  place  there  are  confined  to  the 
material  which  is  already  mixed.  The  flow  rates  and  interconnections  between  the  reactors  forming  the  net- 
work are  determined  directly  from  the  overall  flow  pattern. 

The  finite  difference  phase  of  the  above  algorithm  has  been  applied  to  predict  the  aerodynamic  and  thermal 
patterns  in  a typical  gas  turbine  combustor.  Figures  1,  2,  3,  for  the  case  of  premixed,  prevaporized  fuel  in- 
jection. Solution  of  the  pertinent  governing  equations  yields  the  local  distribution  of  the  velocity  com- 
ponents, p r *e,  turbulence  quantities,  various  species  mass  fractions,  temperature  and  radiation  fluxes. 

Figure  f ’ display  the  evolution  of  the  temperature  field  together  with  the  flux  sums,  defined  as  the 

average  vaii  of  the  two  components  prevailing  at  a point,  by  means  of  longitudinal  and  radial  profiles. 

Particularly  noteworthy  is  the  variation  of  the  x-direction  fluxes  with  maximum  values  lying  close  to  the 
cylindrical  wall.  For  transverse  profiles  a minimum  is  predicted  close  to  the  centre-line,  while  another 
distinct  dip  appears  midway  between  the  combustor  wall  and  the  centre-line.  For  the  particular  radial 
plotting  sections  employed  no  significant  axial  variation  is  detected  in  the  longitudinal  plots. 

The  axial  variation  of  the  y-direction  fluxes  depicted  in  figures  8 and  9 clearly  indicate  the  influence 

of  the  incoming  jets  with  maximum  values  lying  close  to  the  combustor  inlet.  The  transverse  plots  on  the 
other  hand  display  morotonical ly  decreasing  profiles  as  the  wall  is  approached. 

The  z-direction  fluxes  shown  only  for  longitudinal  sections  exhibit  similar  trends  to  the  y-direction 
fluxes . 


Some  Infra-Red  Diagnostic  Techniques  in  Combustion 


Temperature  measurements 

The  measurement  of  surface  temperatures  by  radiant  heat  transfer  is  well  established  in  industry,  for 
example  for  the  measurement  of  glass  and  steel  temperatures,  however  it  is  little  used  in  general  research 
since  it  is  more  expensive  than  a thermocouple  and  depends  on  the  knowledge  of  surface  emissivity.  The  hemi- 
sperical  pyrometer  overcomes  the  latter  problem  provided  the  surface  can  be  accessed  directly. 

Infra-red  imaging  is  a valuable  tool  in  the  design  and  development  of  furnaces  (and  other  heated  structures 
such  as  buildings.)  However,  it  should  be  used  in  conjunction  with  mathematical  modelling  of  the  structure 
heat  conduction  in  order  to  improve  the  precision  of  predicting  the  performance  of  future  designs.  In  parti- 
cular, the  current  swing  towards  ceramic  fibre  insulators  should  be  monitored  carefully  to  check  both  the 
dynamic  response  and  the  effect  of  insulation  anisotropy  on  furnace  design  and  operation. 

With  the  increasing  use  of  coal  in  new  devices  such  as  the  fluidized  bed,  intriguing  problems  are  posed 
concerning  the  respective  temperatures  of  the  gas,  fluidized  solid  inerts,  and  burning  coal  particles.  The 
rapid  response  and  accurate  position  information  of  a two  colour  l.R.  pyrometer  suggest  that  it  could  be  used 
to  resolve  such  problems.  However  current  observations  suggest  that  only  about  10Z  of  the  heat  is  transferred 
by  radiation  in  a fluidized  bed. 


Velocity  measurement 

In  order  to  determine  the  structure  of  two  phase  flames  (i.e.  liquid  or  solid  particulate  fuel),  information 
is  required  on  the  velocity,  cite  and  trajectory  of  the  particles.  In  particular,  the  initial  size  and  velo- 
city of  the  particles  is  required  for  the  mathematical  modelling  of  combustors.  The  gss  velocity  distribution 
is  slso  required  (among  other  parameters)  in  this  hot  hostile  environment,  and  this  can  generally  be  taken  as 
equal  to  the  velocity  of  particles  less  than  1 pm. 

Laser  doppler  anemometers  are  now  used  routinely  to  make  measurements  of  the  velocities  of  particles  in 
flames.  In  this  technique,  the  particles  pass  through  the  intersection  of  two  coherent  beams  of  light  pro- 
duced by  a laser.  In  the  region  of  intersection,  notional  interference  'fringes'  are  formed  which  cause  the 
light  scattered  by  the  particle  to  be  modulated  at  a frequency  proportional  to  their  velocity.  In  addition, 
provided  the  optics  are  suitably  arranged  as  shown  in  Figure  10(5)  the  envelope  of  the  doppler  signal  burst 
can  also  be  used  to  deduce  the  size  of  the  particles.  Although  this  technique  is  normally  applied  with  visible 
light,  there  is  no  fundamental  reason  why  infra-red  radiation  should  not  be  used  similarly. 

An  alternative  technique  for  the  measurement  of  the  size  and  velocity  distribution  of  particles  has  also 
been  developed  at  our  laboratory,  utilizing  a hologram  as  an  intermediate  data  recording  medium ^ . The  holo- 
gram of  the  moving  particles  is  made  using  a double  flash  ruby  laser,  in  which  the  12  n6  flashes  can  be  sep- 
arated by  an  interval  varying  from  1 ps  to  1 ms . Thus  two  holographic  images  of  each  particle  are  recorded, 
separated  by  a small  effective  diatance  depending  on  the  interval  and  the  velocity.  The  interference  fringes 
observed  on  reconstruction.  Figure  11,  can  be  interpreted  directly  in  terma  of  the  velocity  distribution.  It 
is  also  interesting  to  note  that  the  particle  size  distribution  can  be  obtained  from  the  hologram  if  it  is 
placed  in  the  path  of  the  laser  diffraction  particle  size  analyser  previously  developed  at  our  laboratory. 

This  latter  instrument  u6es  s special  photo-detector  and  interfaced  computer  to  rapidly  determine  the  part- 
icle size  distribution  from  the  Fraunhofer  diffraction  of  a collimated  laser  beam  passed  through  the  part- 
icles. Again,  this  normally  uses  a laser  operating  in  the  red  part  of  the  visible  spectrum  732.8  nm  He/Ne. 
However  the  principle  would  hold  for  wavelenghts  increasing  to  the  microwave  band,  with  the  associated  feature 
that  the  upper  size  limit  of  the  device  would  then  extend  to  several  cms . A long  wave  diffraction  device 
such  as  this  could  be  useful  for  continuously  sizing  granulated  materials  as  large  as  several  centimeters. 
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Identification  of  Stirred  Reactor  Networks 

Id  chemical  engineering  plant  decign  procedure,  the  concept  of  the  stirred  tank  reactor  and  plug  flow 
reactor  are  frequently  used  to  identify  reactors  having  strong  mixing  and  feeble  mixing  respectively.  In 
combustors,  the  well  stirred  region  occurs  close  to  the  flame  stabilizer,  and  the  system  decays  to  a plug 
flow  reactor  when  the  turbulence  diasipation  it  complete  a short  distance  from  the  stabilizer.  As  a result 
of  recirculation  and  multiple  jet  entries,  a typical  combustor  can  be  represented  by  a network  of  stirred  and 
plug  flow  reactora  whose  volumes  and  flows  must  be  determined  in  order  to  compute  the  production  of  pollutants 
from  the  chemical  kinetic  rate  equations  for  the  many  individual  reaction  steps  occurring  as  discussed  in  the 
modelling  section  above. 

Infra-red  measurements  can  be  used  to  identify  the  reactor  sequence  by  the  pseudo-random  test  technique 
in  which  disturbances  are  fed  to  the  combustor  then  the  response  at  the  exit  plane  or  various  locations 
within  the  combustor  are  monitored.  Since  only  the  a.c.  response  signals  are  required,  a simple  infra-red 
absorption  measurement  is  adequate  to  monitor  the  response  in  terms  of  concentration  of  CO,  COj , etc.  The 
interpretation  of  the  measurements  is  carried  out  as  follows :- 

The  output  signal  from  a practical  system  is  almost  invariably  corrupted  with  inherent  noise  in  the  system. 
Assuming  the  validity  of  linearity,  application  of  superposition  yields: 

Z(t)  - y(t)  ♦ n(t)  (7) 

where  Z(t),  y(t)  and  n(t)  represent  the  output  signal,  response  signal  and  noise  respectively.  The  mathemati- 
cal expression  relating  the  process  response  and  the  applied  input,  i(t),  is  the  familiar  convolution  inte- 
gral,  i.e.  „ 

y (t ) * ' h(tj)  i (t-tj)  dtj  (8) 

o 

' Here  h(t)  is  the  weighting  function  (or  impulse  response)  of  the  process. 

Combining  equations  7 and  8,  the  measurable  output  to  an  arbitrary  input  signal  is  given  by:  (inputs 
occurring  at  timet  greater  than  T.  in  the  past  are  assumed  to  have  no  effect  on  the  present  output.) 

rT‘ 

Z(t)  • hCtj)  i (t-tj)  dtj  ♦ n(t)  (9) 


The  deconvolution  of  the  above  when  a short  duration  pulse  is  used  yields: 

Z(t)  - Kjh(t)  ♦ n(t)  (10) 

where  K.  is  a constant  determined  by  the  energy  in  the  input  pulse.  The  approximation  becomes  mere  accurate 
as  the  pulse  duration  becomes  shorter.  A step  function  input  yields  the  step  response  which  is  the  integral 
of  the  impulse  response.  The  effect  of  noise  on  the  system  dynamics  represented  by  n(t)  in  the  above 
equation  can  be  substantially  reduced  via  correlation  techniques. 

It  can  be  shown  that  a relationship  exists  between  the  input-output  cross-correlation  function  (ccf) 
and  the  input  auto-correlation  function  (acf),  i.e. 

T 

r 8 

Riz  (t)  * h (tj)  Rjj  (t-tj)  dtj  ♦ Rin  (t)  (11) 
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o 

(c.f.  equation  9) 

where  Rit(i),  Rji(t)  end  Rjn(T)  represent  re.pectively  the  input-output  ccf,  input  acf  and  the  input-noise 
ccf.  Assuming  R— (t ) can  be  represented  by  a delta  function  for  random  inputs,  equation  11  reduces  to  the 
following  siqple  form: 

Riz(T)  ® h (t ) (12) 

The  practical  implementation  of  a random  signal  as  used  in  system  analysis  and  identification  is  a pseudo- 
random binary  signal,  PRBS . . . . , ...  „„„„ 

The  PRBS  approximates  to  a truly  random  signal  due  to  the  properties  of  its  acf.  For  a periodic  PRES,  this 
can  be  shown  to  cosiprise  a series  of  triangular  spikes  of  width  2At  where  At  denotes  the  bit  interval.  Further 
sure,  due  to  its  constant  spectral  density  the  PRES  is  persistently  exciting  - a desired  quality  of  the  input 
signal  for  parameter  estimation  and  identification  studies,  as  the  variance  of  the  estimated  system  parameters 
is  generally  proportional  to  the  inverse  of  the  power  input. 

Amongst  the  generally  available  identification  techniques  including  correlation  methods,  instrumental  vari- 
ables, maximim:  likelihood,  etc,  the  generalised  least  squares,  CIS,  algorithm  has  proven  to  be  statistically 
sufficient  and  is  linear  in  its  formulation.  The  algorithm  can  provide  unbiased  estimates,  a desirable  fea- 
ture, as  data  acquired  for  the  formulation  of  simple  dynamic  models  is  almost  always  corrupted  with  correlated 
noise  and  the  need  to  evaluate  unbiased  estimates  of  the  process  model  and  of  the  noise  becomes  spparent. 
Additionally,  means  exist  to  apply  various  diagnostic  checks. 

The  algorithm  is  applied  to  discrete  input/output  data  sequences  it,  Zt  (fl  N)  to  produce  optimal 

models  of  the  form: 
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where  p it  the  number  of  integrttiont  in  the  tystem  and  noise  models.  W represents  the  forward  shift 
operator  and  k is  the  system  time  delay.  V is  identified  as  the  differencing  operator  (1-W"1)  and  A,B,C  and 
D are  taken  to  be  polynomials  of  the  form: 


A(W_1) 
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B(U_1) 

5 bjW-1  ♦ ... 
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m 
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In  equation  13  the  rational 
extraneous  behaviour. 

transfer 

function 

driven  by  the 

uncorrelated  sequence  £t  comprises 

the  entire 

Equation  13  can  be  re-expressed  upon  removal  of  the  system  integrations  by  data  differencing  in  the 
following  convenient  form: 


AZt  « W'k  Bit  ♦ ~ Ct  (15) 

It  is  worth  noting  that  unless  ADC  * = 1,  the  process  parameter  estimates  will  be  biased.  Hence  the  CLS 
Aims  Co  transform  the  Cera  ADC  iteratively  to  an  uncorrelated  sequence  via  the  following  6teps: 
i)  The  process  parameters  (aj.b^)  are  initially  estimated  by  an  ordinary  least  squares, 
ii)  The  residuals,  viz, 

et  - AZt  " W Bic  (16> 

(>*  denotes  estimated  values) 

are  analysed  to  be  subsequently  transformed  by  autoregression. 
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tii)  The  process  input  and  output  are  filtered  with  the  autoregression  F to  yield: 

F A 
Zj  - FZt 


(17) 

(18) 


iv)  i(  and  Z£  are  modified  by  a new  least  squares  fit  to  comsence  another  iteration  cycle  proceeding  from 
step  ii.  The  final  model  is  expressed  in  the  fora: 


A(FZt)  - W~k  B (Fit)  + Et  (19) 

The  validity  of  the  estimated  process  and  noise  models  given  by  the  CLS  algorithm  have  to  be  established 
by  diagnostic  checks  for  the  statistical  properties  of  the  models.  Currently  available  diagnostic  algorithms 
utilize  model  order  tests  including  determinant-ratio  test,  F-ratio  test,  loss  function  analysis,  pole-zero 
cancellation  and  tests  for  independence.  In  addition,  auto-correlations  of  the  residuals  and  cross-correlaticre 
of  the  input  and  the  residuals  can  be  employed  to  reveal  the  adequacy  of  the  fitted  model. 

Tentative  values  for  the  time  delay  and  model  order  can  be  obtained  via  the  impulse  response  evaluated 
directly  by  cross-correlation  of  the  input  and  output  for  white  inputs,  equation  10.  An  alternative  means 
of  extracting  the  same  information  is  to  use  the  determinant  ratio  tests  The  test  is  especially  valuable 

in  that  it  limits  the  number  of  possible  model  orders  and  associated  time  delays  prior  to  parameter  estimation 

* if  the  signal  to  noise  ratio  is  reasonably  high. 

Loss  function  analysis  and  the  associated  F-ratio  test  are  employed  once  the  process  and  noise  models  have 
been  derived  by  the  CLS  algorithm.  They  are  primarily  concerned  with  the  behaviour  of  the  error  function 
in  the  neighbourhood  of  N0,  taken  to  correspond  to  the  actual  system  order. 

• Additional ly , residual  manipulation  in  terms  of  acf’s  and  ccf’s  provides  a critical  assessment  of  the 
structure  and  order  of  the  estimated  model.  The  a cf  of  the  residuals  is  employed  to  yield  an  indication  of 
the  measure  of  the  whiteness  of  the  residuals.  The  ccf,  on  the  other  hand  establishes  the  statistical  inde- 
pendence of  the  residuals  and  the  process  input. 

The  above  ideas  have  not  yet  been  applied  using  infra-red  diagnostics,  however  using  a tracer  technique 
they  have  been  applied  comprehensively  to  study  the  isothermal  residence  time  distributions  in  the  Lycoming 
combustor  with  the  objective  of  quantifying  the  stirred  and  plug  flow  reactor  locations.  The  predicted 
aerodynamic  patterns  yielded  by  the  finite  difference  procedures  are  utilized  in  conjunction  with  the  CLS 
algorithm  to  determine  the  local  mixing  history.  A detailed  description  of  the  experimental  set  up  and  the 
measuring  environment  is  given  in (10). 

The  investigation  employed  PRBS  and  step  signals  respectively  a6  the  system  mode  of  excitation  for  para- 
meter estimation  purposes.  The  particular  PRBS  adopted  comprises  a period  15,  bit  interval  1.1  sec  and 
amplitude  ♦] volts,  to  a solenoid  valve  introducing  the  tracer.  A number  of  measuring  stations,  represent- 
ative cf  the  distinct  mixing  zones  in  the  combustor,  were  chosen  for  identification  and  parameter  estimation. 

A typical  sequence  of  operational  steps  adopted  in  such  an  analysis  is  suomarized  as  follows: 
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A dttcrminxnt  ratio  teat  waa  employed  to  limit  the  number  of  poeaible  model  orders  for  the  range  of 
time  dclaya  suggested  by  the  deconvoluted  impulse  response.  CLS  parameter  estimation  vas  then  applied  for 
varying  model  orders  and  time  delays.  The  optimum  value  for  the  latter  vas  taken  as  that  which  minimised  the 
sum  of  squared  residuals.  Typically  noise  model  orders  of  10-15  were  adopted  in  an  iterative  scheae  of  5-10 
steps.  The  loss  function  vas  recorded  for  all  the  different  cases  tried. 

Diagnostic  checks  comprising  mainly  the  evaluation  of  the  acf  of  the  residuals,  ccf  between  the  input 
sequence  and  the  residuals  and  pole-zero  cancellation  vere  applied  for  varying  model  orders  and  time  delays. 
Parameter  estimation  and  subsequent  validation  vas  carried  out  in  a strict  iterative  manner.  Finally  the 
impulse  response  function  corresponding  to  the  estimated  model  vas  computed  and  this  vas  compared  with  the 
original  weighting  sequence  obtained  from  the  cross-correlation  analysis  of  the  input/output  data. 

Figure  12  displays  waveforms  labelled  modified  output  and  deterministic  prediction  errors.  These  respect- 
ively refer  to  the  system  response  resulting  from  the  preliminary  data  transformations  (normalising,  trend 
removal,  editing,  etc.)  and  the  computational  mismatch  between  the  predicted  and  measured  outputs. 

The  analysis  reveals  that  the  measuring  station,  immediately  below  the  incoming  primary  jet*,  can  be 
modelled  as  a aecond  order  system  with  a time  delay  of  1.  These  findings  are  quanti tatively  supported  by  the 
excellent  agreement  observed  between  the  predicted  and  measured  pulse  responses.  Figure  13.  Thus  the  system 
response  is  identifiable  by  a transfer  function  having  the  parameters  as: 

»i  - -0.7527  , -0.1179 

b.  - 3.157  , 116.5 

Fire  detection 

It  is  well  known  that  fires  emit  copious  quantities  of  infra-red  radiation,  which,  in  principle  can  be 
used  for  Che  detection  of  unintentional  fires.  A coaroon  problem  with  alarm  systems  of  this  type  is  the 
difficulty  of  distinguishing  between  real  fires  and  false  alarms.  False  alarms  can  arise,  for  exasple  from 
reflected  sun-light,  and  the  electronic  system  associated  with  the  alarm  often  attenpts  to  eliminate  this 
possibility  by  selecting  signals  modulated  at  the  characteristic  flicker  frequency  associated  with  Rayleigh- 
Taylor  instability  of  the  flame.  Nevertheless,  there  are  situations  such  as  sunlight  reflected  from  the 
moving  ocean  surface  where  this  is  not  successful  and  an  alternative  must  be  sought.  The  evolution  of  the 
micro-processor  now  provides  an  important  new  signal  processing  capability  to  detection  systems.  A novel 
application  ia  to  combine  the  signals  from  two  or  more  detectors,  which  scan  the  scene  of  interest  mechani- 
cally, optically  or  electronically,  to  compute  the  two  or  three-dimensional  location  and  size  of  the  fire. 

If  the  'fire'  results  from  s solar  reflection,  or  some  other  source  outside  the  region  of  interest  the  alarm 
is  not  raised.  Such  a system  can  also  be  'taught'  interactively  to  ignore  intentional  fires  or  other  radia- 
tion sources  in  any  particular  section  of  the  overall  region.  The  associated  display  system  can  also  pin- 
point the  location  of  the  fire  and,  if  necessary,  take  selective  action  to  extinguish  it.  The  protection  of 
department  atores,  warehouses  and  even  oil  rigs  by  such  a system  opens  dramatic  new  prospects  in  infra-red 
applications  technology. 

An  interesting  result  from  some  of  our  recent  research  on  oil  tank  farm  fires  shows  that  they  are  strongly 
affected  by  the  air  flow  around  the  tank  resulting  from  cross-wind.  This  is  important  since  the  average 
cross-wind  at  a tank  farm  is  about  6 metres/second.  Since  the  inertial  forces  of  the  air  are  larger  than  the 
flame  buoyancy  force,  the  flame  tends  to  be  drawn  down  into  the  counter-rotating  vortices  in  the  wake  of  the 
tank.  Adjacent  tanks  and  equipment  are  often  located  in  this  region,  and  it  follows  that  they  will  be  sub- 
jected to  ignition  sources  and  strong  convective  heating  in  addition  to  the  radiant  heat  transfer  from  the 
plus*.  As  pointed  out  sbove,  mathematical  modelling  procedures  have  now  advanced  to  the  point  where  these 
complex  three  dimensional  reactive  flow  fields  can  be  predicted,  permitting  accurate  estimates  of  the  hest 
transfer  by  solving  the  relevant  differential  equations  simultaneously  and  integrating  all  the  radiation 
contributions. 

Window  design 

With  the  increasing  use  of  optical  diagnostic  techniques,  we  are  often  faced  with  the  problem  of  window 
design.  One  of  the  simplest  techniques  in  hot  (flames)  environments  is  to  use  an  uncooled  plastic  (plexiglass) 
window.  This  ablates  away  thus  cooling  itself  naturally.  On  a typical  rocket  motor  operating  at  100  ats, 
3000°K  the  window  ablates  about  0.2  sxn/sec.  and  it  is  relatively  simple  and  cheap  to  replace  the  inner  6 mm 
a thick  vindow  after  each  firing. 

If  the  distortion  introduced  by  ablating  windows  cannot  be  tolerated,  e.g.  in  L.D.A.  systems,  then  quartz 
windows  are  recommended.  Often  it  is  necessary  to  protect  the  surface  from  fuel  sprays  and  other  deposits  and 
a screening  layer  of  clean  gas  must  be  used.  In  this  case  it  is  important  that  the  gas  is  introduced  through 
vide , low  speed  jets  at  the  edge  of  the  vindow.  The  jet  will  protect  the  window  for  its  potential  core 
length  ( ' 10  slot  widths)  as  shown  in  Fig.  It*--  Beyond  the  potential  core,  in  the  mixing  region,  the  vindow 
can  easily  become  contaminated  from  the  entrained  material. 

Optical  computers 

With  modern  microcircuits,  the  speed  of  data  processing  has  become  remarkably  fast,  for  example,  current 
state  of  the  art  circuitry  will  permit  a special  purpose  built  circuit  to  obtain  a one-dimensional  Fourier 
transform  at  128  points  in  30  microsecond* . however,  optical  computers  can  operate  at  least  1000  times 
faster  (see  Table  1)  and  for  some  particular  applications  they  offer  distinct  advantages.  The  mathematical 
operations  which  they  can  perform  readily  are  two-dimensional  Fourier  transformation  (as  demonstrated  by  the 
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later  diffraction  particle  tire  distribution  meter  discussed  above),  two-dimensional  cross-correlation  and 
two-dimensional  convolution  functions.  Figure  15  shows  the  typical  arrangement  of  lenses  to  achieve  these 
operations.  Since  the  operations  which  can  be  carried  out  by  an  optical  computer  are  limited,  it  is  usually 
necessary  to  combine  an  optical  computer  with  a digital  processor  to  achieve  an  overall  syatem. 

Combustion  Control 


The  aim  of  an  automatic  mixture  ratio  control  system  is  to  maintain  the  fuel  and  air  supplies  to  a boiler 
or  furnace  in  the  correct  proportions.  The  correct  ratio  between  the  fuel  and  air  is  usually  attained  when 
there  is  just  sufficient  air  to  ensure  complete  oxidation  of  the  fuel  under  the  conditions  existing  in  the 
furnace.  Any  mixture  ratio  control  syatem  should  be  judged  on  the  closeness  with  which  it  approaches  this 
ideal  in  steady  and  unsteady  state  conditions. 

If  mixture  ratio  control  could  be  fitted  to  all  the  combustors  in  the  U.K.,  the  benefit  to  the  country  in 
terms  of  cost  of  fuel  saved  would  be  about  £1000  p.a.  plus  the  advantage  of  mitigating  the  iminent  energy 
shortage.  There  is  therefore  considerable  incentive  to  utilize  or  develop  an  economic  and  reliable  com- 
bustion control  system.  Visits  to  any  but  the  largest  burner  installations  will  soon  show  that  most  present 
monitoring  and  control  instrumentation  is  unreliable  or  unused.  It  is  not  within  the  scope  of  this  paper 
to  survey  the  features  of  all  combustion  control  systems,  however  the  characteristics  of  some  systems  in- 
volving infra-red  devices  will  be  discussed. 

If  the  furnace  is  generally  operated  in  steady  state  the  response  of  the  control  system  does  not  need  to 
be  particularly  fast  and  exhaust  gas  sampling  systems  can  be  used.  Carbon  monoxide  measurements  using  an 
infra-red  absorption  instrument,  usually  of  the  non-dispersive  type,  are  conmonly  used  to  provide  the  control 
signal,  and  the  carbon  monoxide  level  is  maintained  at  some  acceptably  small  value.  The  problems  with  such 
a system  are  largely;  the  need  for  calibration  of  the  instrument;  sampling  problems  due  to  the  fact  that 
carbon  monoxide  concentration  is  far  from  uniform  across  the  flue;  the  problem  of  matching  individual  burners 
on  a multi-burner  installation;  soot,  water  and  other  deposits  in  the  sample  line;  air  in  leakage  in  the 
furnace  or  aa^le  system;  accumulated  soiling  of  the  sample  chamber  windows;  operation  and  'maintainence 1 
by  unskilled  operators;  and  finally,  the  coat  of  proper  maintainence  which  quickly  exceeds  the  cost  of  any 
fuel  saving.  The  overall  experience  and  economics  has  therefore  acted  against  the  installation  of  fuel/air 
ratio  controls  on  smaller  systems  in  spite  of  the  potential  benefit. 

Carbon  dioxide  measurements  can  also  be  based  on  infra-red  absorption,  however  the  concentration  vs. 
fuel/air  ratio  characteristic  becomes  relatively  flat  in  the  vicinity  of  stoichiometric  mixtures  end  it  does 
not  therefore  provide  a very  good  control  signal. 

An  interesting  control  signal  can  be  derived  from  the  flame  radiation.  Since  the  band  radiation  is  pro- 
portional to  the  concentration  of  the  emitting  species  multiplied  by  approximately  the  fourth  power  of  the 
temperature,  both  of  which  peak  close  to  atoichiometric  mixtures,  one  has  a signal  with  a pronounced  peak  as 
a function  of  fuel/air  ratio,  and  a peak-seeking  feedback  control  system  can  easily  locate  this  peak  even 
in  the  presence  of  noise.  In  the  case  of  a natural  gas  flame  the  major  emitters  are  COj  and  HjO  whose  infra- 
red bands  carry  most  of  the  radiated  energy.  Provided  the  mixing  time  in  the  flame  is  small  compared  with 
the  kinetic  time,  which  is  true  for  all  high  intensity  gas  burners  for  both  pre-mixing  or  nozzle-mixing 
systems,  the  peak  is  quite  reliable.  However,  if  the  flame  is  operating  as  a slow,  non-premixed  diffusion 
flame,  then,aa  with  a candle  flame,  the  air  fuel  ratio  cannot  be  defined,  and  a peak  seeking  system  cannot 
be  used. 

Kith  oil  or  coal  flames,  most  of  the  radiant  energy  is  due  to  luminous  soot  in  the  flame  and  the  radiation 
depends  primarily  on  the  fourth  power  of  the  soot  (i.e.  flame)  temperature  and  non-linearly  on  soot  concen- 
tration. In  this  case,  again  neglecting  diffusion  type  flames,  the  peak  in  radiation  occurs,  however  it  is 
now  at  the  soot  emission  boundary.  That  is,  if  the  flame  is  operated  any  leaner  the  flame  temperature  drops 
due  to  excess  air,  whilst  any  richer  results  in  the  flame  temperature  dropping  due  to  the  emission  of  un- 
burned  agglomerated  soot  from  the  burner.  The  radiation  peak  therefore  usually  corresponds  to  the  desirable 
operating  point  of  the  burner.  This  point  is  usually  closer  to  stoichiometric  when  the  burner  is  operating 
at  full  load  than  when  it  is  operating  at  low  load. 

It  is  interesting  to  note  that  a peak  seeking  detection  system  does  not  need  calibration,  is  insensitive  to 
dirty  windows,  and  can  readily  identify  individual  burners  in  a multi-burner  installation. 

With  the  advent  of  cheap  digital  processors,  it  is  now  possible  to  operate  the  burner  at  any  required 
air/fuel  ratio  away  from  the  peak,  with  only  occasional  excursions  to  the  peak  to  obtain  its  location.  Such 
systems  can  also  attend  to  light-up  and  operate  several  system  simultaneously.  An  interesting  festure  is  thst 
they  can  'remember*  and  periodically  update  (as  necessary)  the  optimum  air  and  fuel  valve  settings  for  any 
.demanded  heat  output  and  thus  an  optimally  fast  feedforeward/feedback  control  characteristic  can  be  obtained 
for  dynamic  load  following. 

Conclusions 

b'e  feel  a burning  enthusiasm  for  digital  processors  in  combustion  diagnostics  and  foresee  that  the  future 
is  particularly  bright  for  infra-red  systems. 
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fig.  1.  Predicted  Radial  Velocity  Profiles 

(Axial  Cross-section)  for  the  Lycoming 
Combustor.  (Isothermal,  K « 4) 
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Fig.  2.  Predicted  Axial  Velocity  Profiles 

(Axial  Cross-Section)  for  the  Lycoming 
Combustor,  (Isothermal,  K * 4) 


Fig.  3.  Predicted  Radial  and  Tangential  Profiles 
(Transverse  Cross-Section)  for  the 
Lycoming  Combustor,  (Isothermal,  I ■ 13) 
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Fig.  12.  Estimated  Model  for  Position  13 

Together  with  Auxilliary  Waveforms. 
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Fig.  13.  Comparison  Between  the  Predicted 
and  Measured  Pulse  Responses  for 
Position  13. 
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